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Abstract 

An important tool for the analysis of results of numerical simulations of lat¬ 
tice QCD is chiral perturbation theory. In Wilson chiral perturbation theory 
the effects of the hnite lattice spacing a are taken into account. In recent 
years the effects of isospin splitting on the masses of hadrons have been in¬ 
vestigated in Monte Carlo simulations. Correspondingly, in this article we 
derive the expansions of the masses of the pseudoscalar mesons in chiral per¬ 
turbation theory at next-to-leading order for twisted mass lattice QCD with 
three light quark flavours, taking the mass difference between the up and 
down quarks into account. The results include terms up to orders in the 
quark masses, Arrd in the mass splitting between up- and down quarks, and 
a? in the lattice spacing, respectively. 
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1 Introduction 


For investigations of QCD by means of numerical simulations on a lattice, chiral per¬ 
turbation theory has become an indispensable tool. It allows to describe the dependence 
of hadron masses and matrix elements on the quark masses in terms of expansions in 
powers of the quark masses, modified by logarithmic terms. On the one hand this of¬ 
fers the possibility to extrapolate the results of numerical simulations with larger quark 
masses into the physical regime with small up- and down-quark masses. On the other 
hand, comparing the quark mass dependence observed in Monte Carlo calculations with 
the formulae from chiral perturbation theory yields the values of fundamental low-energy 
constants of QCD, the Gasser-Leutwyler coefficients. For reviews see DEI- 

The discretisation of space-time on a lattice with lattice spacing a implies a breaking 
of chiral symmetry and leads to discretisation effects in observable quantities. In the 
case of lattice QCD with Wilson fermions the latter are generally of order a. A way to 
reduce lattice artifacts is to employ the formulation of twisted mass lattice QCD, which 
implements a chirally rotated mass term ism]. Setting the twist angle to the value 7i/2 
(“full twist”) implies automatic 0{a) Symanzik-improvement [HE]- The ETM collab¬ 
oration has done simulations with two mass-degenerate quark flavours, representing the 
up and down quarks, and with 2 -|- 1 -|-1 quark flavours, meaning a degenerate doublet of 
up and down quarks and a non-degenerate doublet of charm and strange quarks, with 
full twist in both sectors EEIE]. 

Chiral perturbation theory for twisted mass lattice QCD has been developed in m 
for N{ = 2 quark flavours. It has been extended to W = 1 + 1 non-degenerate 
quarks and N{ = 2 -|- 1, including the strange quark, in [H], and to Af = 2 -|- 1 -|- 1, 
with a degenerate doublet of up and down quarks, in na, and in US] in the form of 
a charmless chiral Lagrangian. The resulting formulae for meson masses reveal the 
dependence on the lattice spacing a in the form of additional terms proportional to 
powers of a. Chiral perturbation theory taking lattice effects for Wilson fermions into 
account is called Wilson chiral perturbation theory. 

In recent years the isospin splitting produced by the small mass difference between 
up and down quarks has come into the focus of numerical investigations of lattice QCD 
[I71IIB1II2I. This splitting contributes to the mass splittings between charged and neutral 
hadrons. An additional source for these mass splittings are electromagnetic effects, which 
have been studied in the framework of lattice QCD -|- QED in [TSl [2ni ED CH] • 

The purpose of this paper is to take the mass splitting between the light quark masses 
into account in Wilson chiral perturbation theory for twisted mass lattice QCD with 
Af = 1-|-1-1-1 quark flavours. We calculate the masses of the members of the pseudoscalar 
meson octet to next-to-leading order (NLO). The small isospin mass splitting Am = 
— TTiu is included to order Am, and lattice effects are included to order a^. 

The purpose of this article is twofold. On the one hand, the results include the limit 
of a vanishing twist angle, which is ordinary untwisted lattice QCD with three non¬ 
degenerate Wilson quarks. For this case our calculations reveal the effects of the lattice 
discretisation on the isospin splitting for pseudoscalar meson masses. In addition to 
that, even the further limit of vanishing lattice spacing a is of interest, as it includes 
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the effects of mixing in next-to-leading order, which are not completely included in the 
original work [22] , 

On the other hand, simulations of twisted mass lattice QCD including strange quarks, 
are in practice done by adding a fourth heavy charm quark and implementing full twist 
in both the up-down and the strange-charm sector in order to ensure automatic 0{a) 
improvement in all observables. For the case of the masses of pions, which are made of 
up quarks and down quarks only, it turns out that full twist in the up-down sector alone 
is not sufficient for 0(a) improvement, and that the fully twisted strange-charm sector is 
necessary. Our calculations make the influence of unaccompanied strange quarks on the 
pion masses explicit and in this way points to the need of including the charm quarks 
in the simulations. 


2 Chiral perturbation theory 


In twisted mass lattice QCD the twist in the quark sector has to be implemented ortho¬ 
gonally to the mass splitting in order to maintain positivity of the fermion determinant 
[a Eg. In this letter we introduce the mass twist in the u-d-sector via the third Gell- 
Mann-matrix A3. Consequently the mass splitting is implemented by means of Ai. The 
strange quark is introduced as a standard Wilson fermion. The three-flavour quark mass 
matrix is thus given by 


M = ml + {ms — m) 


/1 

1 \ 


/ m 

4 Am 


(-1 - 

1 

-I- -AmAi = 

4Am 

m 

0 

V3 

73 V 

2 

0 

0 

ms) 


^ . m^ + md . 

m = -, Am = mj — m^. 


The chiral rotation entailing the twisted mass is given by 
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/m -h i75/i 

4 Am 
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M exp 
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Mexp 
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m - iys/i 







1 0 

0 

ms/ 


H = m sin(ci;), m = m cos(a;), u = arctan 



( 1 ) 

( 2 ) 

(3) 

(4) 


The effective Lagrangian describing the physics of the pseudoscalar mesons is con¬ 
structed via the power-counting scheme of [21], accounting for terms of order {m, a} in 
leading order (LO), and {m^,ma, a^} in next-to-leading order (NLO), respectively. The 
Lagrangian is formulated in terms of the matrix-valued held 


U = exp 


L^O a=l 


(5) 


It contains the meson helds (pa and the low-energy constant (LEC) Fq, which equals the 
pion decay constant to leading order. The quark masses and the lattice spacing a enter 
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p = Pol = 2Woal , 


the Lagrangian through the variables 


X = 2BoM , 


( 6 ) 


where Bq and Wq are additional LECs. To leading order the chiral effective Lagrangian 
reads 

C^o = ^ Tr {d,U d,U^) - ^ Tr + Ux^)-^Tt (p + Up^) . (7) 

In numerical simulations of twisted mass lattice QCD a chirally rotated mass matrix 
like Eq. ([3]) has to be used. On the other hand, in chiral perturbation theory it is more 
convenient to employ the physical basis. It is obtained by applying the inverse of the 
above chiral rotation, such that the mass matrix retains its original form Eq. ([T]), and 
the lattice term goes over to 


p{u) = exp 


, OJ 



/P + iP3 

0 

nl 

1 A 3 — 

L 2j 

pexp 

■ 1 

lAs — 

L 2j 

0 0 

p-ip3 

0 

Po/ 


= ^ (2p + Po) 1 - ^ (Po - P) As + ipsAs , 


p = Po cosfcnj, 


Ps = Po smfcuj 


( 8 ) 

(9) 


In order to set up chiral perturbation theory for the computation of meson masses and 
other observables, the effective Lagrangian has to be expanded around its minimum. In 
contrast to the case of continuous space-time, where the minimum is located at vanishing 
meson helds, = 0, its position is shifted for twisted mass lattice QCD and is located 
at 


= 0 for i ^ 3, 
and the corresponding value of the held U is 


03 — 


FqPs 
X + P ’ 


Uq = exp 



( 10 ) 


( 11 ) 


It has turned out that a suitable parametrisation of the held U, taking the expansion 
point Uo into account, is given by 

U = (12) 


such that the minimum of the Lagrangian is now at V 
shifted meson helds 


U' = exp 


0 L 
[ 1^0 


1. Expanding in terms of the 


(13) 


with the help of the Baker-Campbell-Hausdorh series Eza has the advantage that spuri¬ 
ous three-point vertices are eliminated and that just one quartic vertex remains in the 
Feynman rules, which simplihes the loop calculations a lot. 
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In order to derive the meson masses in NLO, the effective Lagrangian has to be 
extended by terms of the next higher chiral dimension. Neglecting external fields and 
constant terms, this leads to the NLO lattice yPT Lagrangian [21l[ID] 




- Li [Tr {d^Ud^U^)] - L2 [Tr [d^UdM^)] - L3 Tr [[d^Ud^U^]^) 

+ L4 Tr Tr (x^f/ + U^x) + ^4 Tr Tr (p^U + U^p) 

+ L5 Tr (d^[/d^[/HUx^ + xf/^)) + W5 Tr (d^Ud^UHUp^ + 

- Le [Tr (x'^t/ + U^x)f - We Tr (x^U + U^x) Tr (p^U + U^p) 

- L, [Tr (x^f/ - f/tx)]' - IT7 Tr (x^l 7 - U^x) Tr {p^U - U^p) 

- Ls Tr (x^Ux^U + U^xU^x) - Tr (^p^Ux^U + U^pU^x) 

- [Tr (ptf/ + U^p)f - it; [Tr (p^t/ - p)^ 

- Ws Tr (^p^Up^U + U^pU^p) . ( 14 ) 


The masses in NLO receive tree-level contribntions from £nlo and one-loop contribu¬ 
tions from £lo- The loop contributions produce divergences that require renormalisation 
of the low energy coefficients Li, Wi and IT/. 

The necessary steps for deriving the meson masses are as follows. The mass matrix ([T]) 
and the twisted lattice term ([S]) are inserted in the Lagrangian ([TT|) . and the resulting 
expression is expanded in the meson fields (pa up to linear terms in order to find the 
minimum. The minimum of the Lagrangian at NLO is given by 


I ^ F2 + 8IT6(2x + Xs) + 8IT8X 

0P3 -HfiPTO-v -U -vh-v 4- 1 firo(''v 2 4 . 


-U n(r,^\ 


( 15 ) 


As this quantity enters the calculation only quadratically, its terms ~ a? can be neg¬ 
lected. Next the parametrisation of the fields according to Eq. fll2p is implemented, such 
that the minimum of the Lagrangian occurs at the origin of the shifted meson fields. 
Expanding around this minimum up to quadratic terms in the fields gives us the inverse 
meson propagator at tree level. To LO the inverse propagator reads 


iG ‘(p^) =pn + OTlo . 


( 16 ) 


The LO mass squared matrix DUlo contains non-diagonal terms that generate mixings 
in the Ti^-p sector and among the kaon fields. The loop calculations require the LO 
propagator and the quartic vertices. The vertices are obtained by expanding the NLO 
Lagrangian up to quartic terms in the shifted meson fields, which leads to complicated 
expressions containing momentum dependent terms. 

At NLO the inverse meson propagator receives additional tree-level contributions and 
self-energy terms from the loop diagrams, 

iG~^{p‘^) =p'^l + Tli^o-A + P^B . ( 17 ) 
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The matrices A and B are rather complicated and contain various off-diagonal mixing 
terms. The divergent contributions in the self energies can be separated from the finite 
ones and absorbed into the NLO renormalised LECs LJ via a modified minimal sub¬ 
traction scheme as shown in [22]. Field renormalisation can be achieved by writing the 
inverse propagator as 


iG-^(p^)= 1 + 


B 


p^l + OHlo — a — -{B, OJTlo} 


1 + 


B 




(18) 


such that the expression in square brackets contains the normalised kinetic term. The 
squared meson masses are then obtained by diagonalising the matrix [OJIlo — A — 

i{s,ajtLo}l. 


3 Results for the isospin splittings 


To LO in the limit Am 0 the method outlined above reproduces the known results 

m 


= 2BQm -j- 2aWo cos (a;) -|- 


a^Wo^sin^(uj) 


Bom 


+ 0(a^ 


m]^ = Bo(m + ms) + alTo(l + cos(a;)) -|- 


la^W^sin^icj) 


™?=3 


2Bo{m -\- 2ms) + 2alTo(2 -l- cos(a;)) -|- 




Bofh 


+ 0{a-^ 


(19) 

( 20 ) 
( 21 ) 


We list these expressions here since they are used as convenient abbreviations in the 
NLO results. 

With our choice of the generators for the quark mass splitting and twist, the neutral 
pion is TTi = 7 r°, and the charged pions are linear combinations of the mass eigenstates 
712 and TTs. The charged pion 7 r 2 represents the component that does not mix with the 
7] and is “orthogonal” to twist (A 3 ) and to mass splitting (Ai), and is therefore least 
affected by these factors. The mass of the charged pion 712 is given in NLO by 
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=^l + 


32Bl 




32a2iy2 


{-2LI- LI + 4Ll + 2Ll)7h^ + ( - + 2Ll)mm, 

32aBoWo + ( - ^4 + wi)m, cos(a;) 

+ [-2LI- LI- 2WJ - 1^5^ + 4W^ + 2W^)7ficos{uj) 

- 2W^ - + LWq + 2W^" +[-Wl + 2 Wq) cos(a;) 

+ {-ALl - 2Ll + 4W^ + 2W^ - AWq - 2W'^") sin^{uj) 


— 9 

mi 


32n^F^ 


LI - 2Ll + 2W^ 




\ 

r rUgl 

• 9 / \ 



sm^(a;) 

/ 

[2rfiJ 



A 2 


A 2 


+ 0{Am\a^^ 


( 22 ) 


In NLO the pions appear to lose their automatic 0(a)-improvement at full-twist, due 
to the presence of the term proportional to {—Wl + hhg)a. This term stems from the 
presence of the strange quark and is independent of u. 

We also note that in the pion masses dependencies on the strange quark mass rrts 
show up, that are independent of u and persist in the continuum limit. As the pions 
are made from up and down quarks, this term proportional to mg is counterintuitive. 
In the context of QCD in the continuum it has its origin in the fermion determinant or, 
in terms of Feynman diagrams, in strange quark loops, that lead to NLO terms in the 
effective Lagrangian. The associated LECs L\ and Lg have very small empirical values, 
however [2S]- 

In the calculation of the NLO pion masses for non-degenerate quark masses we find 
that the mass splitting between the neutral pion and the charged pion 112 is induced by 
the TT^-r] mixing, and surprisingly appears to be at least of order Am?. These terms are 
manageable in the tree-level contributions, but lead to quite complicated dependences in 
the loop contributions, in particular in the 0{a^) terms. Since the quark mass difference 
Am^ is very small, we decide to list the pion mass splittings including the Am^ terms up 
to order a only, neglecting 0(a‘^) terms. Then the mass difference between the neutral 
pion TTi and the charged pion 712 is 
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^ alTo(l — coscu) 
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.' [(-4 - 3Wi 

F^{m^ - m) L 

+ {LI + LI + 20U^ + 4Ll + 2Wl + W; - 3W^ - lOW^ - 


-3L\ - 2Ll + QLl) 


Fi \ms 

3Wl - + 3W^) 

2 Wg){l — cosu) 

{3LI + LI + ALl + - 3W^ - 2W^){1 - cosu) 



1927r2Fo2 0 




- 0{Am^, a^). 


(23) 


We observe that at full twist, cos a; = 0, also terms of order a are present. They are not 
proportional to mg. The presence of these terms is a consequence of the fact that the 
strange quark is not a member of a twisted doublet. Namely, if the strange quark would 
be accompanied by a fourth quark, forming together a fully twisted doublet, automatic 
0 {a) improvement would forbid the existence of such terms. 

Due to the mass-twist the charged pious represented by the helds 112 and also gain 
different masses, which leads to the splitting 


Ml - Ml = -^ 

773 772 


0 


{-4LI - 2Ll + 4Wg + 2Wg - 4Wi' - 2Wg) sin^)^) + 0(0?). 

(24) 

This mass difference is a pure lattice artifact that vanishes in the continuum limit a —?■ 0. 

We have also calculated the mass splitting between the kaon masses in NLO including 
terms of order a^. As the characteristic features show up at order a already, we display 
the mass splitting between neutral kaons Kq, Kj and charged kaons it's, including 
terms linear in a, which is 

























— M‘^± = Bo Ami 1 


16^0 
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^ — 2ij4 — Lg + 4 ijg + ^U^rh + ^ — Z74 — Z75 + 2ijg + 2 L\ 
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3271^F^ 
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-^In 

3 

1 
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Bo{ms — m) 2 
C>(Am^, a^). 


mlln -fhl\n(^ 


A 2 


A 2 


1 - 


aWo{l — coso;) 
Bo{ms - m) 


(25) 


It vanishes in the limit Am = 0, as it should. In the continuum limit the expression 
reduces to that given in 122 ]. 

In contrast to the case of pions, the kaon mass splitting has contributions linear in 
Am, part of which is produced by pion and eta loops. It is thus signihcantly more 
sensitive to the tiny quark mass difference Am than the pion mass splitting. Therefore 
the kaon masses represent observables ideally suited for numerical investigations of the 
quark mass difference Am in lattice QCD. For studies with three flavours of Wilson 
quarks the formula shows the lattice artifacts proportional to oAm. 

We have also obtained the expression for the mass of the r] meson in the case of 
non-degenerate quark masses. It is quite lengthy, and as it does not reveal any new 
characteristic features compared to the pion masses, we refrain from presenting it here. 


4 Conclusions 

To summarise, the masses of the pseudoscalar mesons for the case of non-degenerate up, 
down and strange quarks have been obtained in chiral perturbation theory for twisted 
mass lattice QCD including lattice terms quadratic in the lattice spacing a and in Am. 
Automatic 0{a) improvement of the pion masses at maximal twist in the up-down sector 
is lost. At maximal twist, terms of order a are present, which have two different origins. 
On the one hand, at NLO the pion masses get an equal contribution of order a that stems 
from the presence of the strange quark in the effective Lagrangian. On the other hand, 
the inverse propagator contains off-diagonal terms that represent a mixing between the 
neutral pion and the rj meson. These yield a contribution to the neutral pion mass which 
is proportional to Am?. Apart from continuum terms (a°) it has terms of order a that 
spoil automatic 0{a) improvement at maximal twist. We present the mass splitting 
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between neutral and charged pious including terms of orders Am^ and a. In addition, 
the twisted mass generates a mass difference between the charged pions of order a^. 

Finally, the lattice corrections to the mass splitting between neutral and charged 
kaons, which are proportional to Am, are presented including terms of order a. In the 
kaon sector there is also mixing, in this case between and Kq and between and 
Kj. As for the pions, it leads to terms of order a that spoil automatic 0(a) improvement. 

If the strange quark were introduced with a partner (charm quark), associated with 
another twist in this heavy sector, the lattice artefacts of order a and higher would 
look quite different. In particular, full twist in both sectors would ensure absence of all 
0(a) terms. Our results make the idea that such inclusion of a twisted fourth quark is 
advantageous for reducing unwanted lattice artifacts explicit. 
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